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Abstract 

In 2006, the fourth author of this paper proposed a graph-theoretic model of interface 
dynamics called competitive erosion. Each vertex of the graph is occupied by a 
particle that can be either red or blue. New red and blue particles alternately get 
emitted from their respective bases and perform random walk. On encountering a 
particle of the opposite color they kill it and occupy its position. We prove that on the 
cylinder graph (the product of a path and a cycle) an interface spontaneously forms 
between red and blue and is maintained in a predictable position with high probability. 



Figure 1: Competitive erosion spontaneously forms an interface. In the initial state S{0) 
each vertex of the 30 x 30 cylinder graph is independently colored blue with probability 
.33 and red otherwise. After 325 time steps of competitive erosion, an interface has formed 
between red and blue. 
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1 An Interface In Equilibrium 


We introduce a graph-theoretic model of a random interface maintained in equilibrium by 
equal and opposing forces on each side of the interface. Our model can also be started from 
a heterogeneous state with no interface, in which case an interface forms spontaneously. 
Here are the data underlying our model, which we call competitive erosion: 

• a finite connected graph on vertex set V ; 

• probability measures fii and fi 2 on V ; and 

• an integer 0 < fc < — 1. 

Competitive erosion is a discrete-time Markov chain {S{t))t>o on the space of all subsets of 
V of size k. One time step is defined by 

S{t + l) = {S{t)U{Xt})-{Yt} (1) 

where Xt is the first site in S{tY visited by a simple random walk whose starting point 
has distribution ^i, and Yt is the first site in S{t) U {Xt} visited by an independent simple 
random walk whose starting point has distribution p, 2 - 

For a concrete metaphor, one can imagine that S{t) and its complement represent the 
territories of two competing species. We will call the vertices in S{t) “blue” and those 
in S{tY “red”. According to (1), a blue individual beginning at a random vertex with 
distribution p,i wanders until it encounters a red vertex Xt and inhabits it, evicting the 
former inhabitant. The latter returns to an independent random vertex with distribution 
p ,2 and from there wanders until it encounters a blue vertex Yt and inhabits it, evicting the 
former inhabitant. 

If the distributions fj-i and /i 2 have well-separated supports, one expects that these dy¬ 
namics resolve the graph into coherent red and blue territories separated by an interface 
that, although microscopically rough, appears in a macroscopically predictable position with 
high probability. The purpose of this article is to prove this assertion in the case of the 
cylinder graph Cyl„ = Cn x Pn, where Cn is a cycle of length n and Pn is a path of length 
re. Here x denotes the Cartesian (box) product of graphs. We identify Cn with (iZ)/Z 
and Pn with (^Z) n [0,1]. We take pi and p 2 to be the uniform distributions on Cn x {0} 
and Cn X {!}; respectively: blue particles are released from the base of the cylinder and red 
particles from the top. 

Taking k = an^ for some 0 < a < 1, it is natural to guess that the stationary distribution 
of the Markov chain S{t) assigns high probability to the event 

A,^n ■■= [S C Cyl„ : x [0, a - e] C 5 C x [0, a + e]}, (2) 

that is, the event that sites below the line y = a — e are all blue and that sites above the 
line y = a + e are all red. Our main result confirms this guess. 

Theorem 1. Given e > 0 there exists a positive constant d = fi(e) such that 

vrn(A,n) > 1 - 6“"^" (3) 

where Tin is the stationary distribution of the competitive erosion chain on Cyl„. 
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The proof of Theorem 1 will show that the predicted interface forms quickly, even if the 
initial state 5(0) is heterogeneous (such as a checkerboard of red and blue, or a random 
initial state like the one shown in Figure 1). 

Theorem 2. Let = inf{t : S{t) £ Ae,n}- For any e > 0 there exist positive constants 
c = c(e), d = d{e), N = N{e) such that for all n > N and all subsets S C Cyl„ of cardinality 
[an^J we have 

P(r, > dn^ I 5(0) = 5) < 

See the last section for some related models exhibiting interface formation. We also 
mention a superficially similar process called “oil and water” [4] which has an entirely 
different behavior: the two species do not form a macroscopic interface at all. 

1.1 Comparison with IDLA 

Internal diffusion limited aggregation (IDLA) is a fundamental model of a random interface 
moving in a monotone (outward) fashion. IDLA involves only one species with an ever¬ 
growing territory 

Iit+l)=Iit)U{Xt} 

where Xt is the first site in I{ty visited by a simple random walk whose starting point 
has distribution /ri. Competitive erosion can be viewed as a symmetrized version of IDLA: 
whereas I{t) and I{t)^ play asymmetric roles, 5(t) and S{t)'^ play symmetric roles in (1). 

IDLA on a finite graph is only defined up to the finite time t when I{t) is the entire vertex 
set. For this reason, the IDLA is usually studied on an infinite graph and the theorems 
about IDLA are limit theorems: asymptotic shape [11], order of the fluctuations [1,2,6], and 
distributional limit of the fluctuations [7]. In contrast, competitive erosion on a finite graph 
is defined for all times, so it is natural to ask about its stationary distribution. To appreciate 
the difference in character between IDLA and competitive erosion, note that the stationary 
distribution of the latter assigns tiny but positive probability to configurations that look 
nothing like the predicted horizontal interface. Competitive erosion will occasionally form 
these exceptional conhgurations: for example, at a tiny but positive fraction of times t the 
boundary of the set S{t) is a vertical line! The proof of Theorem 1 is delicate since there 
are so many exceptional configurations: in terms of cardinality, the desired set Ae,n is an 
exponentialy small fraction of the set of all recurrent configurations. 

1.2 Idea of the proof 

To prove that the stationary distribution concentrates on the small set A = Ae,n, we 
identify a Lyapunov function h on the state space which attains its global maximum in A 
and increases in expectation 

E(/i(5(l)) -/i(5(0))) > a > 0 (4) 

provided 5(0) is sufficiently far from A. The function h is as simple as one could hope for: 
the sum of the heights of the red vertices. The heart of the proof is Theorem 5, which uses an 
electrical resistance argument to establish the drift (4) for a suitable notion of “sufficiently 
far from A”. Since the function h{a) is bounded by n^, we then use Azuma’s inequality to 
argue that the process h{S{t)) spends nearly all its time in a neighborhood of its maximum. 
This establishes Theorem 4 (a statistical version of Theorem 1). 
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The main remaining difficulty lies in showing that if ^(O) is sufficiently close to A, then 
the chain S{t) hits A quickly with high probability, establishing Theorem 2. This is done 
using stochastic domination arguments involving IDLA on the cylinder. These ingredients 
together with a general estimate relating hitting times to stationary distributions (Lemma 4) 
establish Theorem 1. 

1.3 The level set heuristic 

Before restricting to the cylinder we mention a heuristic that predicts the location of the 
competitive erosion interface for well-separated measures , /i2 on a general finite connected 
graph, which we assume for simplicity to be r-regular. Let g he & function on the vertices 
satisfying 

= Iti - Ai2 (5) 

where A denotes the Laplacian 

^9^^) '■= \ ^i9{x) - g{y)) 

and the sum is over vertices y neighboring x. Since the graph is assumed connected, the 
kernel of A is one-dimensional consisting of the constant functions, so that equation (5) 
determines g up to an additive constant. 

The boundary of a set of vertices S' C F is the set 

dS = {v £ V — S : u ~ tc for some ui G S}. (6) 

Consider a partition of the vertex set C = Si U B U S2 where dSi = B = 882- (Think of 
Si and S2 as the blue and red territories respectively, of the sort we might expect to see in 
equilibrium, and the sites in their common boundary B have indeterminate color.) 

Let gi for i = 1, 2 be the Green function for random walk started according to gi and 
stopped on exiting Si. These functions satisfy 


Agi = Hi on Si, 
gi = 0 on Sf. 

The probability that simple random walk started according to Hi first exits S, at x G B is 
-A5j(x). 

To maintain equilibrium in competitive erosion, we seek a partition such that Agi ~ Ag 2 
on B, that is 

A(gi - 52) ~ /Ui - M 2 - 

(Exact equality holds except on B.) Thus by (5), the function g — {gi — 52 ) is approximately 
constant. Since gi vanishes on B, the equilibrium interface B should have the property that 

g is approximately constant on B. 

The partition that comes closest to achieving this goal takes 5i to be the level set 

= {x : g{x) < K} (7) 
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for a cutoff K chosen to make = k. An application of the maximum principle shows 
that for this choice of Si, the maximum and minimum values of — {gi — g' 2 ) differ by at 
most 

max \g{x)-g{y)\, 

x£Si,yfSi,X'^y 

suggesting that the right notion of “well-separated” measures yi and ^2 is that the resulting 
function g has small gradient. 

1.3.1 Mutually annihilating fluids 

The divisible sandpile [13,14] is a deterministic analogue of IDLA. In a competitive version 
of the divisible sandpile, the edges of our graph form a network of pipes containing a red fluid 
and a blue fluid. The two fluids are injected at respective rates gi and ^2 and annihilate 
upon contact. One can convert the above heuristic into a proof that this deterministic model 
has its interface exactly at a level line of g (where g is interpolated linearly on edges). 

1.4 Diffusive sorting 

We define here a process called diffusive sorting, introduced by the fourth author, which 
couples the competitive erosion chains on a given graph for all values of k, using a single 
random walk to drive them all. We will not use the coupling in this paper, but find it to 
be of intrinsic interest. 

It is convenient to imagine that the blue and red random walks start from vertices and 
V 2 respectively that are external to the finite connected graph G and have directed weighted 
edges into G: each edge {vi,v) has weight gi{v). 

The state space of diffusive sorting consists of bijective label lings A of the vertex set of G 
by the integers {1,..., iV}, where N is the number of vertices. We imagine that the labels are 
detachable from the vertices and can be carried temporarily by the random walker. Consider 
an initial labeling Aq. A random walker starts at vi carrying the label 0. Whenever it comes 
to a vertex whose label exceeds the label the walker is currently carrying, the walker swaps 
labels with that vertex, dropping its old (smaller) label and stealing the new (larger) label 
for itself. Eventually the walker will visit the vertex labeled N and acquire the label N for 
itself, at which point we may stop the walk, since the particle will carry the label N forever 
after and no labels will change. The vertex labels are now {0, ...,N — 1}. We now increase 
all those labels by 1, and call the result Ai/ 2 - (We have noticed that this process bears a 
strong resemblance to an algorithm in algebraic combinatorics called promotion [16], but 
we have not explored this connection.) To get from A 1/2 to Ai, we apply the same process 
from the other side, releasing a random walker from V 2 bearing the label A^-|- 1, and letting 
it swap labels with any vertex it encounters whose label is smaller than its own, until it 
acquires the label 1; then we decrease the labels of all vertices by 1, obtaining Ai. 

Diffusive sorting is the Markov chain {Xt)t>o on labelings, each of whose transitions from 
Xt to Xt+i is as described in the previous paragraph. To see how this chain relates to 
competitive erosion, for each integer fc = 0,..., — 1 let 

Sk{t) = {v : Xtiv) < k}. 

Then it is not hard to check that has the law of the competitive erosion chain. (The key 
observation is that the vertices labeled 1 through k in Aq must be a subset of the vertices 
labeled 1 through A; -|- 1 in A 1 / 2 , and similarly for going from A 1/2 to ^i-) 
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The following result shows that the stationary distribution of diffusive sorting on the 
cylinder graph is concentrated on labelings that are uniformly close to the function n^y. 


Theorem 3. For each e > 0 there is a constant d = d{e) > 0 such that 


A : sup 

\{x,y) 
r.2 y 

>'] 

[ ix,y)£Cy\„ 


J 


^—dn 


where nn is the stationary distribution of the diffusive sorting chain {Xt)t>o on Cyl„ 


The level set heuristic (§1-3) makes a prediction for diffusive sorting on a general graph: 
if the gradient of g is not too large and we label the vertices by {g{v) : v £ V} instead of 
by {1,..., A^}, then 7r„ should concentrate on labelings not too far from the function g. 

We mention, as an aside, that if instead of alternating between blue walkers and red 
walkers we use walkers of just one color, we get a sorting version of ordinary IDLA. For any 
k > 1 and t > 0 the law of Sk{t + k) is that of the IDLA cluster with k particles. (This is 
trivial for k = 1 and the rest follows by induction on /c.) 


2 Connectivity properties of competitive erosion dynamics 

Implicit in the statement of Theorem 1 is the claim that competitive erosion has a unique 
stationary distribution. In this section we formally define the competitive erosion chain and 
prove this claim. 


2.1 Formal definition of competitive erosion 

The graph Cyl„ is a discretization of the cylinder (M/Z) x [0,1] with mesh size To avoid 
degeneracies, we will always assume n > 2. Let Pn = H [0,1] be the path graph with 
edges ^ ~ for A: = 0,1,..., n — 1. Let Cn = be the cycle graph obtained by 

gluing together the endpoints 0 and 1 oi Pn- Let 


Cyl„ = CnX Pn 


( 8 ) 


with edges (x, y) ~ (x', y') if either x = x' and y ~ y', or x ~ x' and y = y'. We also add 
a self-loop at each point in Cn x {0} and Cn x {!}, so that Cyl„ is a 4-regular graph. We 
will use Cyl„ to denote both this graph and its set of vertices. 

For t = 0,1, 2,... let (Ai*^)s>o and be independent simple random walks in 

Cyl„ with 


P(aW = (x, 0)) = - = = (x, 1)) 


h+5) 


n 


for all X G Cn- That is, each walk starts uniformly on Cn x {0} and each walk 
starts uniformly on x {!}. Given the state S{f) of the competitive erosion chain at time 
t, we build the next state S{t + 1) in two steps as follows. 


S(*+2) = S{t)U{W,} 


where T(t) — inf{s > 0 : xi** ^ S{t)}. Let 


s(t+i)=s(t+i)-{F;+4(} 
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where T{t + 5 ) = inf{s > 0 : G S{t + ^)}. 

To highlight the symmetrical roles played by S{t) and its complement, we will often think 
of the states as 2 -colorings rather than sets: let 


Also define for i = 1 , 2 


crt{x) 


1 X e S{t), 

2 X ^ S{t). 


Bi := Bi{a) := {x G Cyl„ : a(x) = i}. 


(9) 


( 10 ) 


Thus 

S{t) = Bi{at). 

We will use color 1 and “blue” interchangeably and similarly color 2 and “red”. Now for 
each a (the proportion of blue sites) strictly between 0 and 1 , competitive erosion defines a 
Markov chain on the space 

n := {cr G {l,2}*^y^",#{x G Cyl„ : a{x) = 1} = L«|CylnlJ}- (11) 

We also set 

n' := {a G {1, 2}^^^",#{x G Cyl„ ; a{x) = 1} = L«|Cyl„|J + 1}. (12) 

A single time step of competitive erosion consists of a step from to 0,' followed by a step 
from n' back to 

2.2 Notational conventions 

We will often use the same letter (generally C, D, c or d) for a constant whose value may 
change from line to line (or sometimes even within one line); this convention obviates the 
need for distracting subscripts and should cause no confusion. 

For any process and a subset A of the corresponding state space r(A) will denote the 
hitting time of that set. When a; is a state in the space and A is a set of states, we will 

write {t{A) < K \ uj} to denote the event that, starting from cj, the process hits A in at 

most K steps. Finally, in all the notation the dependence on n will be often suppressed. 

2.3 Blocking sets and transient states 

Definition 1. Call a subset A C Cyl„ blocking if Cyl„ \ A is disconnected and the subsets 
Cn X {0} \ A and Cn x {1} \ A lie in different components. 

We next prove that the competitive erosion chain has exactly one irreducible class and hence 
has a well-defined stationary measure. We start with a definition. 

Definition 2. For two disjoint blocking subsets A,B C Cyl^j we say that A is over B if 

1. A and Cn x {0} \ B lie in different components of Cyl„ \ B, and 

2 . B and Cn x { 1 } \ A lie in different components of Cyl„ \ A. 
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Cn X {1} 



Cn X {0} 

i. a. 


Figure 2: i. The configuration cj*. ii. A. configuration with a blue blocking subset over a 
red blocking subset. 

Lemma 1. The competitive erosion chain has exactly one irreducible class. Moreover any 
a € it that has a blue blocking set over a red blocking set is transient. 

Proof. Consider the configuration cr* in Figure 2 where the lowest [an(n + 1)J vertices are 
colored blue. To prove the first statement notice that from any a one can reach u*. Since in 
the target configuration there is exactly one blue component, we look at the closest vertex 
with a valne 2 from Cn x {0}. Since there is a blue path from Cn x {0} to that point the 
Markov chain allows us to change it to 1 and similarly at the other end. Thus we are done 
by repeating this. 

To prove the second statement we prove that starting from (Tq one cannot reach any 
configuration like the one on the right that has one blue blocking set over another red 
blocking set. We formally prove this by contradiction. Let u be a configuration with a bine 
blocking set B over a red blocking set A. Now the competitive erosion chain evolves as 

<70) 0-1/2,0-1,0-3/2 • • • 

where for every non-negative integer k, ak € ^ and o-fc_|_i /2 G Assume now that there is 
a path 

o'o, 0-1/2, 0-1,0-3/2 ..., (Tt = 0-. 

Let t' and t" be the last times along the path such that at least one vertex of A is bine 
and similarly B contains at least one red vertex respectively. If such a time does not exist 
let ns call it — 00 . Since o-q does not have a blue blocking set over a red blocking set, 

max(r',r'') > — 00 . 

t' must be a half integer (since at half integers there is one more blue particle) and similarly 
t" must be an integer. Thus 

/ / // 

T ^ T . 

Next we see that we cannot have r' < r". For, suppose otherwise. Then for all times greater 
than t', A \s entirely red. Thus no blue walker crosses A at any time greater than t' . So i? 
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must already be blue at t' and stays blue through till t implying that t" < t' . By a similar 
argument we cannot have t' > t" . Hence we arrive at a contradiction. □ 

2.4 Organization of the proofs 

In Section 3 we state a weak “statistical” version of the main result, Theorem 4, and provide 
its proof using hitting time estimates whose proof appear in Section 5 . A key step in the 
proof, the construction of a suitable Lyapunov function, appears in Section 4. This section 
uses the theory of electrical networks. A short review of a few basic facts about electrical 
networks that are assumed in this section appears in Appendix B. In Section 6 we deduce 
the stronger Theorem 1. As a simple corollary we obtain Theorem 3. This last section 
uses an estimate for IDLA on the cylinder proved in Appendix A. We also apply a hitting 
time estimate for submartingales several times throughout the paper. The proof of the 
estimate is included in Appendix C. The final section discusses some related models and 
future directions. 

3 Statistical version of the main theorem 

Given e > 0 we define the set 

:= {a e Q : |(T“^(l)A{y < q:}| < en^} . (13) 

where A denotes the symmetric difference of sets, AaB := (A^ n B) U (A n B^). 

Theorem 4. Given e > 0 there exists c > 0 such that for large enough n 


TTniQe) > 1-6 


This theorem says that, after running the competitive erosion chain for a long time, one 
is unlikely to see many blue sites below the line y = a oi many red sites above it, which 
(see Remark 1 hi. below) is a weaker statement than Theorem 1. 

3.1 The height function 

We now introduce a function on the state space that will appear throughout the rest of the 
article and will be used as a Lyapunov function, along lines sketched in (4). The interested 
reader can learn more about Lyapunov functions in [5]. Given a G (recall 11 and 12), 

let us define the height function of a as 

h{^)= C^-y) ( 14 ) 

{x,y)eBi{a) 

where Bi{a) is defined in (10). Glearly 

sup h{a) = h{(Tmax) = «(! - + 0{n), (15) 

where amax is any element of a with a{x,y) = 1 for the lowest [an{n + 1)J vertices. Also 
notice that for all t 

\h{at) - h{at-i)\ <1. (16) 
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Figure 3: Examples of the various good sets defined in Table 1. Fg is defined for technical 
purposes and is not included in the figure since it essentially looks like by Remark 1. 


Set 

Description 

Defined in 


all sites outside an e-band are the right color 

(2) 

Ge 

at most en^ sites are the wrong color 

(13) 

F, 

height function within en^ of its maximum 

(17) 


at most en rows have sites of the wrong color 

(25) 


Table 1: Four different kinds of good sets. 


Definition 3. Throughout the rest of the article given u G U we say that a vertex 
V = {x, y) has the wrong color if either y < a and a{v) is red or if y > a and a{v) is blue. 

Definition 4. For e > 0 let 

F^ := {cj G O : h{a) > (a(l - - e)n^}- (17) 

where h{-) is the height function defined in (14). 

Remark 1. One easily observes the following inclusions: 

1. r,2 cGecr^ 

2. A C A 

3. A C Q2e- 

The following lemma states that starting from any configuration a the process takes 
O(n^) steps to hit F^. 

Lemma 2. Given any small enough e > 0 there exist positive constants c = c{e),d = 
d{e), N = N{e) such that for all n > N and a G Q 

> 2dn^) < e-^^\ ( 18 ) 
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The next lemma asserts that once the process has hit Tg it is likely to stay in a neigh¬ 
borhood T 2 e for a long time. 

Lemma 3. Given e > 0 there exist positive constants c = c{e),d = d{e),N = N{e) such 
that for all n > N and a G T^ 

P,(r(Q \ Tse) > e"”') > 1 - e-^^\ (19) 

The remainder of this section deduces Theorem 4 from Lemmas 2 and 3. The Lemmas 
themselves are proved in Section 5. 


3.2 Bounding stationary measure in terms of hitting times 

We now prove a general lemma on Markov chains relating hitting times to the stationary 
distribution, which will be used in the proofs of both Theorems 1 and 4. Roughly, it says 
that if one can identify subsets A C B oi the state space of a Markov chain such that (1) 
A is hit quickly from all starting points, and (2) it takes a long time to escape B from all 
starting points inside A, then the stationary distribution assigns high probability to B. This 
along with Lemmas 2 and 3 would then allow us to conclude that r 2 e has high stationary 
measure. 


Lemma 4. Let u}{-) be an irreducible Markov ehain on a finite state space Ai. Suppose 
Ac B C Ai. Let ti,t 2 ,Pi,P 2 be sueh that 


max{Pt^(r(^) > ti)} < pi 

min{P^(r(i?‘^) > t2)} > 1 - P2- 

loGA ^ 


Then 


tt{B^) < 


— +Pl +P2, 

t2 


where vr is the stationary distribution of the Markov chain u}{-) on Ai. 


( 20 ) 

( 21 ) 


Proof. First we note that by the ergodic theorem for Markov chains for all cu G Al 




lim 

t^OO 


C B-) 
t 


( 22 ) 


The proof of the above is standard. For example, it directly follows from [12, Theorem 4.9]. 
We now claim that the following bound is true: 


C 

sup E l{ui{k) G B’^) <ti + {pi +P2){t2)- 

^ 6 -^ k=l 


To see this decompose the sum on the left hand side: 


*2 t{A)M 2 t2 

E^^l{co{k)cB^)=E^ l{co{k) C B^)+E^ Y l{oj{k) C B^). 

k=l k=l k=T{A)/\t2 
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The first sum is at most ti on the event that t{A) < ti and at most t 2 on the complement. 
The second sum is at most t 2 on the event that after t{A) the process exits B in fewer than 
t 2 steps and 0 otherwise. Hence adding these up we get the upper bound 

t2 

sup ^ l{uj{k) G B'^) <ti+ pit 2 + P 2 t 2 - (23) 

We now decompose the sum £ B^) into blocks of length t 2 , i.e. 

t t2 2t2 

k=l k=l k=t2-\-l 

The theorem now follows from (22) by using (23) for each of the above sums on the right 
hand side. □ 


3.3 Proof of Theorem 4 

Notice that it suffices to prove that for given small enough e for all large enough n 

7r(r2e) > 1 - (24) 


for some c = c(e) > 0. By the lower containment in Remark 1 this proves that 


The proof now follows immediately from Lemma 4 with the following choices of the param¬ 
eters: 


A 

B 

h 

h 

Pi 

P2 


r. 

Ls, 

div? 

^d2n 

e-cin^ 

g-C2n2 


where ci,C 2 ,di,d 2 are chosen such that the hypotheses (20) and (21) 
satisfied. Lemmas 2 and 3 allow us to do that. 


of Lemma 4 are 

□ 


4 The main technical result 

We now build towards the proofs of Lemmas 2 and 3. The next result is one of the key 
technical ingredients of the paper. Recall (11) and (12). Given u G H U for any a G 
{0, , 1} let us say that the line y = a is a bad level if there exists x = {x, a) such that 

(t(x) has the wrong color. We now define the set of configurations a that have very few bad 
levels. 
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Figure 4: The sets Ri and R 2 of Definition 6 . Note that the two blue islands and the red 
island are not included in either set. 


Definition 5. Given e > 0 for all n let be the set of all conhgurations a € Q such that 

#{y : o-((Ty)) = 2 for some (j,y) G Cyl„, with y < a} < en and (25) 
#{y : cr((j, y)) = 1 for some (j, y) G Cyl„, with y > a} < en. (26) 

Define D), similarly. 

Recall the competitive erosion chain at defined in (9). This section is devoted entirely to 
the proof of the following theorem: 

Theorem 5. (Lyapunov Function) Given e > 0 there exists a = a(e) > 0 such that for all 
n sufficiently large and for all u G 

E(/i(cJi) — h{ao) \ (Xq = cf) > a 

where h{-) is the height function defined in (14). 

For A C Cyl„ and v G Cyl„ dehne 

i/A(^^) = E^[y^(^c)] (27) 

where , yr{A‘=) ) is the point at which the random walk exits A and is the expectation 

under the random walk measure on Cyl„ with starting point v. 

Definition 6. Given cr G D U D' define Ri = Ri{a) C Gyl„ to be the set of all points of 
color 1 in Gyl^ reachable by a monochromatic path of color 1 in a from a point in Cn x {0}. 

Similarly let R 2 = R 2 {x) be the set of all points of color 2 in Cyl„ reachable by a 
monochromatic path of color 2 from a point in Cn x { 1 }. See Figure 4. 


Lemma 5. For u G D, 


E{h{ai) - h{ao) \ ao = a) 


1 

n 


Ecto 


uiec„x{i} 


^ X] ^RiDo)(^) ( 28 ) 

1JeCnX{0} 

(29) 
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where 


<^0) <7i/ 2, o-i 

form the steps of the competitive erosion chain. 

Proof. The first walker of color 1 starts uniformly on Cn x {0} and stops on first hitting a 
site of color 2, that is, on first exiting iii(cr). Writing 

/i(cri) - h{ao) = [h{ai/ 2 ) - h{ao)] + [h{ai) - h{ai/ 2 )], 

we see that the expected value of the first term in brackets on the right side is 

Dec„x{o} 

Now a walker of color 2 starts uniformly on Cn x {1} and stops on first exiting 
so the expected value of the second term in brackets is 



n 


CO 


w£CnX{l} 


- 1 


(this term appears with an expectation because ai /2 is random). 


□ 


4.1 Energy and flows 

Next we relate the expression on the the right side of (28) to energy of flows. We first recall 
some standard notation in electrical network theory. On a graph G = (E, E) let u) ~ u 
signify that rc is a neighbor of v. Also let E be the set of directed edges where each edge in 
E corresponds to two directed edges in E, one in each direction (except for self-loops, which 
correspond to just one directed edge in ill). A flow is an antisymmetric function / : iH —)• M 
(i.e. a function satisfying f{w,v) = —f{v,w)). Note that by definition the value of a flow 
on a self-loop is 0. Define the energy of a flow / by 

= \ Z (30) 

{v,w)gE 

For any flow / : iH —>■ M define the divergence div / : F —>■ M by 

dWf{v) = '^f{w,v). (31) 

W'^V 

Note that for any flow 

^ div / = f{x, y) + f{y, x) = 0 

XGV x,y£V 

y^x 

since / is antisymmetric. For disjoint subsets A,BcV and a flow / we say that the flow 
is from A to ii if div f{x) = 0 for all vertices x except vertices of A and B and the sum of 
the divergences across vertices of A and B are non-negative and non-positive respectively. 
For more about flows see [12, Chapter 9]. 

Now we state a key lemma relating Lemma 5 to energy of flows on Cyl„. 
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Lemma 6. For any A C Cyl„ such that A n {Cn x {1}) = 0, 

i ^ HA{v) = ini£{f) 

Dec„x{o} ■' 

where the infimum is taken over all flows from A P| {Cn x {0}) to A^ such that for {x, y) & A 

div(/)(i,!,) = | » 

Similarly for any B C Cyl„ such that B n {Cn x {0}) = 0, 

E HB{v) = if£if) 

where the infimum is taken over all flows from {Cn x {1}) CiB to B'^ such that for {x, y) £ B 

div(/)(i,!,) = | » 

The proof is deferred to Appendix B. We now briefly sketch how it helps in proving 
Theorem 5. 

By (28) and Lemma 6 one sees that 

E(/i(c 7 i) - h{ao) I fjo = cr) = 1 - £{fi) - E[£{f 2 )], 

where fi and /2 are flows on Ri{a) and -R 2 (c’'i/ 2 ) satisfying the divergence conditions men¬ 
tioned in Lemma 6 and have minimal energy. Note that the second term on the RHS has 
an extra expectation to average over the random intermediate random configuration cri/ 2 - 
It turns out that the flow with minimal energy on Ri is the “natural” flow induced by 
random walk started uniformly on Cn x { 0 } and killed on exiting Ri and similarly for i? 2 - 
Theorem 5 is now proved by estimating £{fi) and £{f 2 )- Consider the voltage function on 
Cyl„ that is linear in the height. £{fi) can be interpreted as the voltage difference between 
Cn X {0} and the boundary of Ri with glued boundary condition. The above heuristics are 
formalized in Appendix B. Now if the boundary is not horizontal enough, gluing vertices 
at different potentials causes a voltage drop. We estimate this voltage drop en route to 
proving Theorem 5 by constructing suitable flows on Ri and i? 2 - 

4.2 Towards the proof of Theorem 5 

For technical purposes we introduce a few definitions. Recall Ri and R 2 from Definition 6 . 

Definition 7. Given a £ D U D' let Ri be the set of all points in Cyl„ reachable by a 
monochromatic path of color 1 in cj from a point in Cn x { 0 } which does not hit the set 
Cn X {!}. 

Similarly define R 2 . See Figure 6 . 

For i = 1,2 define 

n = T{Rfi) (32) 

to be the exit time from Ri. Also let 

fi = t{R\) = Ti a T{Cn X {!}) (33) 

and similarly f 2 . 
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Figure 5: The green and yellow segments depict and y” respectively for k £ Cn- 


Remark 2. Recalling the notation in (27) we see that 

HR^(a)iv) < 

since on the event fi < ti, y^^ = 1 and hence 


Vri < Vfi = 1 . 

For notational simplicity we write 

Hl(0 = «!,.(•) = 

ftO) = «2,4-) = 

We use the first notation when the underlying a is clear from context. 
For every k £ Cn, a G fl U 12' let 

yfc(o') = inf{y : {k,y) ^ Ri} 
yr(<7) = sup{y : (fc,y) ^ ^ 2 }- 


(34) 

(35) 


See Figure 5. 

As usual we will suppress the dependence on a and denote Vk Vk* 

respectively when a is clear from the context. 


Lemma 7. For cr G 12 U 12' 


n 


Y.n,(k,o) < It. 


k^Cn 


Vk 


k^Cn 


n 


Tn2(m > It yi 


k^Cri 


kSCn 


(36) 

(37) 


Proof. We will only prove (36) since the arguments for (37) are symmetric. By Remark 2 
it suffices to prove that 

^i(A:,0)< yl (38) 

k^Cn k^Cfi 
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i. a. Hi. 

Figure 6; i. A general configuration with Ri and R 2 . Note than in ii. the shaded red region 
and the top light-blue line (i?i nCn, x {!}) are outside Ri. Similarly the shaded blue region 
and the bottom orange line are outside R 2 . Notice the difference between Ri and Ri. In 
in. d°^^Ri is represented by the purple curve. The rest of dRi (bottom left) is represented 
by the red curve. 

Note that Ri satisfies the hypothesis of Lemma 6. Thus using the lemma we will prove the 
above by finding a flow f* from Ri n {Cn x {0}) to R^ such that for all {k, y) £ Ri 

div{f*){k,y) = I ” ^ 

^ I 0 otherwise 

and 

Sin 

Let 

f*i{k,y),ik,y+-)) = --foi {k,y) £ Ri, 0 < y < yl (40) 

n n 

and let every other edge have zero flow. It is easy to check that f* satisfies the required 
divergence conditions and (39). See Figure 7 i. □ 

The next lemma can be thought of as a weaker version of Theorem 5. It shows that at 
all times the height function has an “almost” non-negative drift whereas Theorem 5 states 
that outside fie the drift becomes significantly stronger. 

Lemma 8. For all a £ 0, 

E(/i(cti) - h{ao) I do = cr) > -- 

n 

where h{-) is the height function defined in (14). 

Note that the above lemma is sharp for do = dmax (15), up to a constant factor. 


E 

n n 


Vk- 


(39) 
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Proof. By Lemma 5 


E{h{ai) - h{ao) \ ao 


a) = -E 
n 


WSCn ^ {1} 


I Y 

Dec„x{o} 


Thus by Lemma 7 


E(/i(cri) - h{ao) I do 


a) > —E 
n 


Y yV^^i/2) 

kSCn 


I Y 

/c£Cn 


( 41 ) 


where yl{-) and yl*{-) are defined in (34) and (35) respectively. We make the following two 
easy observations; For all do = d G fl and k € Cn, 

yn^i/2) > yn^o) (42) 

yr(c'o)-yfc(do) > (43) 

n 

The first one is a straightforward consequence of the fact that di/2 has one more blue vertex 
than dQ. To see the second one, notice that for all k G 4z/reZ, either y^(do) = 0 or 

{{k,0),{k,yl{ao) - i)} G .Ri. 

In the first case (43) follows since yl*{cro) > 0. In the second case, Ri and i?2 are disjoint 
and therefore yl.{crQ) — 4 ^ R 2 . Hence by definition (43) follows. 

Thus 


E(/i(di) — /i(do) I do = d) > —E 

n 


Y yti^i/ 2 ) 


k^Cn 


IY yfc(^o) 


k^Cn 


> - 


Y yk*i^o) 


kSCn 

1 


kSCn 


> 

n 


(44) 

(45) 

(46) 


The first inequality is (41), the second and third follow immediately from the above two 
observations. Hence we are done. □ 


4.3 Bending flows 

We prove Theorem 5 in this subsection by improving upon Lemma 8 which was proved 
using the trivial flow on the cylinder graph that sends all the mass along the vertical edges. 
However if the boundary of Ri (Definition 6) is not horizontal then flows that send some 
mass along the horizontal edges are more optimal. This is because the horizontal boundary 
edges have significant harmonic measure and the random walk started from Cn x {0} has a 
significant chance of exiting Ri through these edges. Thus to prove Theorem 5 we modify 
the trivial flow by sending some mass along the horizontal edges. See Figure 7. We begin 
with a few definitions. 
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Figure 7: Bending the trivial flow to improve upon Lemma 8. In i. every vertical edge carries 
mass In ii. the flow carries mass ^ up to height yk where it splits into ^ (vertically) and 
1 — ^ (horizontally). The dk’s are chosen to be ■ Also note the flows are split only 

for the columns where and y” are close to each other to beat (43) used in Lemma 8. 

For any m G (0,1) let us define the following sets: 

W = {£:yr > m} (47) 

IF' = {£:ye < m}. (48) 

Also let us define the following properties for c > 0 and m G (0,1): 

Property 1. There exist cn/A horizontal lines below the level y = m — c/4 that intersect 
Cyln\Ri. 

Property 2. There exist cn/A horizontal lines above the level y = m + cjA that intersect 

Cyl.\i?2. 

The statement of the next lemma roughly says that if there are sufficiently many bad 
levels below or above a certain horizontal line then the height function experiences a strong 
drift. The proof uses the idea of bending flows outlined above. 

Lemma 9. Given c,d>0 suppose a € 0, is such that one of the following holds: for some 
m, 

i. Property 1 holds for c and \ W\ > dn. 
ii. Property 2 holds for c and |IF'| > dn. 

Then for large enough n, 

E(/i(cJi) — h((To) \ ctq = a) > a (49) 

where a = a(c, d). 

Proof. Let m be as in the statement of the lemma. Because of the obvious symmetry 
between the two conditions we will only discuss the proof of case i. Without loss of generality 
we can assume c = d. 
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y = m 

y = ye 



Vi) 


Figure 8: Proof of Lemma 9 in case of Property 1. 


Let Wi be an arbitrary subset of W of size cnj^. We choose an arbitrary pairing of points 
in W\ with horizontal lines mentioned in Property 1, i.e. for every I in W\ we associate a 
distinct horizontal line y = yi that intersects dRi and yi <m — cj Also let 

^k-,ye) (50) 

be the point on the line y = yi closest to the point [l^yi] such that {jt^yi) G dRi. exists 
since the line y = ye intersects dRi (see Figure 8). Now for all i G Wi, 

(51) 

(53) 

c 1 
8 n 

The extra — - term follows from (43). 

Also for all i G IL2, 

y*e-ye>^- (54) 

We now consider the two following sub-cases: 


as VFi C W. 
Let 


Hence for all ^ G Wi \ W2 


Thus we have 


yT -yi> c/4 

W2 = {£ G Wi : %* > m - 
c 

Ve -ye> g- 
1 


n 




** - vt] > - 
n 


k^Cn 


E 


k&Wi\W2 


(i) IIF2I < \Wi\/2. 


Then since |bFi| = ^ and hence |bFi \ TT2I > by (53) we have 


1 

n 


E isr 

k^Cn 



1 

n 


(55) 
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Thus in this case by (44) 


1 

E(/i(cri) - h{ao) I CTO = a) > —- 

96 n 


and hence we are done. 


(ii) |W^2| > \Wi\/2. 


The proof in this case is more involved than the previous case. We first claim that in this 
case there exists a flow fi from Cn x {0} H Ri to Cyl„ \ Ri such that for all (/c, y) G Ri 


div(/i)(/c,?/) = 


^ ify = 0, 


0 otherwise 


and 


£(/0 S - E 


* 

Vk 


kSCn 


1536 


(56) 


(57) 


Before proving the above claim we show why it suffices and implies (49). By taking the 
set A in Lemma 6 to be Ri followed by using Remark 2 we get 

1 
n 


Hl,ao{v) < Y 


(58) 


wec„x{o} 


1 

< - 
n 


1JeCnX{0} 

yl 

k^Cn 


1536 


By Lemma 5 


E(h(cJi) - h{ao) \ uq = a) = ^E[ Y 7f2,^i/2('w)] - ^ Y (59) 

1 
n 


The first term is bounded from below by — vTi^o) (as in the proof of (44)) and the 

kec 

second term by (58) is upper bounded by 

1 yU<^o) 


kSCn 


1536 


Hence we get 

E(/i(cJi) - h{ao) I do = ct) > 




> 


k^Cji 

A 1 


kSCn 


+ 


1536 


1536 n 

Note that the term in the brackets is greater than — ^ by (43). Thus the proof of Lemma 
9 is complete in this case except for the proof of the initial claim. We now prove the claim 
by defining the flow /i as shown in Figure 8. Recall (52). For k ^ IT2 let 

fi{{k,y),{k,y + i)) = for all {k,y) G Ri and y < yl- 
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For k G W 2 


fi{{k,y),{k,y +-)) = 

n 


-i if 


if yk<y< vl 


y <yk 


/i(0',yfc), O' + -,yfc)) = --——\lj = k,k + -,...jk 

n n n 

where jk was defined in (50) and 0 < < 1 are to be specified later. /i is 0 on all other 

edges. It is easy to see that the flow fi satisfies the divergence conditions (56). 

To see why (57) is true first observe that 


S{fi) < - ^ [yk + d\{y*k — yk) + (1 — dfc)^] + — ^ 


Vk- 


k&W2 k^W2 

This is because by construction: 

• for k ^ IT 2 the flow along the line a: = A: is ^ for ny'^ edges. 

• for A G IT 2 

— the flow along the line x = A: is ^ for nyk edges and ^ for n{y’l, — yk) edges. 


— the number of edges with flow 


1 dh 


on the line y = yk is at most n. 


Plugging in dk = we get that the expression on the RHS equals 


E 

fceW 2 

E 

k&W2 


yk + 


{y*k - Vk) 


1 


{y*k - yk) + ij n 
yl + (yk — yV) + 


+ -Y.y*k 

k0V2 
iy*k - Vk) 


(yfc-yfc) + ij n 


- 1 . yk’ 


k^W2 


- Y. yk + +- Y. y*k + - Yl 

k&W2 k0V2 kGW2 

{yl - Vk? 


{yk - yl) + 


{yl - yk) 

{yl - Vk) + 1 . 


1 

n 


T.yi-iE 


k^Cn 


keW2 


.{yl-vk) + 1_ 


Now since |VFi| > ^ and we are considering the case IIT 2 I > 2 ^) have IIT 2 I > 

Also for k G IT 2 by (54) 


cn 
12 • 


g < 2/fc - 2/fc < 1- 


Plugging these in the expression 


■ (y’k-Vkf , 


E 


gives us (57). Thus the proof of Lemma 9 for case i. is complete. 


□ 
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The proof of Theorem 5 involves considering a few cases and showing that the hypotheses 
of Lemma 9 are satisfied in each case. We start with some definitions. 

Let Cyl^ be the slightly larger graph Cn x Pn +2 where we identify the path Pn +2 with 

{- : o = -1,0,1 ,... ,n + 1}. 
n 

Recall the definition (6) of the boundary of a set of vertices. The following is a technical 
definition needed for convenience (see Figure 6): 

Definition 8. Let dRi be the boundary of Ri U [Cn x {^}) in Cyl((. Define d°^^Ri to 
be the subset of dRi that is visible from Cn x {1}, i.e., the set of all points in dRi that 
are connected to x {1} in the subgraph induced by Cyl„ \Ri- As usual we have similar 
definitions for (9i?2, d^*‘R 2 - 

Note that since Ri n {Cn x {1}) = 0, 

dRi C Cyl„. 

Define the graph Cyl* to be the graph Cyl„ along with the additional diagonal edges 

(\n nJ \ n n /J 

n n J \ n n J 

for alH = 0,... n — 1, J = 0 ... n — 1 where the addition in the first coordinate is in /Z. 
Call a subset B C Cyl„, ^-connected if it is a connected set in the graph Cyl*. 

Remark 3. ft follows from [17, Lemma 2] that for i = 1,2, d°^^Ri is a ^-connected set. 

Let 

y = {y {k,y) G for some k G Cn}- (60) 

Note that T is a connected subset of [0, ^ ..., 1] since is ^-connected. Also for every 

k G Cn let 

y{k) = mi{y : {k,y) G d°^^Ri} (61) 

y^'^^ = sup{y:(fc,2/)G5°"*Ri}. 

See Figure 9. Recall the definition of y^ and y” from (34),(35). Note that for every k G Cn, 

yl<y{k)- (62) 

We make another simple observation: For any k G Cn there exists i G {k — k,k + ^} such 
that 

y(D_i<y«. (63) 

n 

This follows since the point {k,y^^'^) G d°^^Ri and hence either = 0, or (/c,y^^)) has a 
neighbor colored 1 in Cyl„. 

We are now ready to prove Theorem 5. As we remarked earlier it suffices to show that 
the hypotheses of Lemma 9 are satisfied. 

Proof of Theorem 5. By hypothesis cr G D \ Also observe that for any ct G D one of 
the following three conditions holds: 
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Figure 9: The yellow segment depicts for k € Cn- The green segment depicts for 
those k for which ^ y{k)- Note that unless d°'^^Ri intersects the column 

X = k ai more than one point. Also we include the purple segment at the bottom and 
the black segment at the top depicting and y|^ respectively to illustrate the differences 
between the four definitions. 

(i) supT — inf T > e/ 20 . 

(ii) supT — inf T < e /20 and | supT — a\> e /10 
(hi) supT — inf T < e /20 and | supT — ct| < e /10 

where the set y was defined in (60). 

Case (i): Let m = ^(sup y + inf y),c = ^. By hypothesis sup y — inf y > c. We first show 
under this hypothesis both Properties 1 and 2 hold. 

By Remark 3, y is connected. Thus there exist cn/4 values in y less than m — cjA. In 
other words: 


Property 1 holds with m and c as defined above. 

Again using connectedness of T, cn/4 values in y are greater than m + c/4. Also notice 
that for every k < supT the line y = k must intersect Ri since the line y = supT has a 
vertex lying on the boundary of Ri and hence we have a monochromatic path of color 1 
from the level y = 0 to y = supT— Thus there exist cn/4 horizontal lines above the level 
y = m + c/4 that intersect Cyl„ \ i? 2 - Hence 

Property 2 holds with m and c as defined above. 

To verify the rest of the hypotheses of Lemma 9 we have to show that either the set W or 
W' is large. Recall y^^^ from (61). Clearly at least one of the inequalities 

\{k : y*-^^ > m}| > n/2, (64) 

|{A: : y^^^ < m}\ > nj^ (65) 

is true. Recall from (47) and (48) 

IT = {£ : yf > m} 

W = {i-.y\< m} 
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When (64) holds, by (63) 


Similarly when (65) holds by (62), 

IW'1>1 

Thus all the hypotheses of Lemma 9 are satisfied in case (i) and hence we are done. 

To verify the hypotheses of Lemma 9 in the remaining cases we start by making an 
observation. Let 

W” := >miy--} 

n 

W" := <supT}. 

Now we have 

l^"l > ”, (66) 

\W"'\ = n (67) 

(the first follows from (63) and the second from (62)). Thus by choosing m = inf y or supT 
the conditions for sets W and W' respectively in the hypotheses of Lemma 9 are always 
satished . 

The remainder of the proof shows that in cases (ii) and (iii) either Property 1 or 2 holds 
by choosing m = inf y or supT- This would then complete the proof by Lemma 9. 

Both cases (ii) and (iii) have two subcases, arguments for which are symmetric. We will 
discuss only one of the cases. In case (ii) we will discuss the case 

supT — a> e/10 


and skip the case supT — a < —e/10. In the case supT — cr > e/10 


inf T — a > —-— 

“ 10 20 


e 


Now this implies that there are at least red vertices below the line y 


inf T — ^ since 


inf y 


40 ^“ + 


e 

40' 


In particular there exist at least ^ disjoint horizontal lines below the line y = inf y — ^ 
with at least one red vertex, i.e. Property 1 is true for m = infT and c = ^. 

In case (iii) we use the fact that ct G O \ fig and hence there are at least en horizontal 
lines with a red vertex below the line y = a or en horizontal lines with a blue vertex above 
the line y = a. We will discuss only the first case. Now by hypothesis in (iii) 


inf T > a 



Hence there exist at least ^ disjoint horizontal lines below the line y = inf T — ^ with at 
least one red vertex, i.e. Property 1 is true for m = inf y and c = ^. 

The symmetric arguments which we skip will show that Property 2 holds in the other 
cases with m = supT and c = ^. Thus we are done. □ 
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5 Proofs of hitting time results 


In this section we provide the proofs of Lemmas 2 and 3 . We first state the following general 
lemma about hitting times of submartingales. The statement involves a few parameters and 
can be slightly difficult to follow. However it will be useful in subsequent applications. Let 
u}{t) be a stochastic process taking values in an abstract set D. Also let 5 : 2 ? —)■ M be a 
real-valued function. Let Tt be the hltration generated by the process uj{t) up to time t. 
Also define Xt := g{uj{t)). 

Lemma 10 . Let Ai, ^2,01 > 0 . Suppose 

|g'(a;)| < Ai for all oj ^ V 

\Xt — Xt-i\ < A2forallt. 

Also suppose that B <ZT) is such that for any time t 

E{Xt - Xt-i\Xt-i) > ( 70 ) 

Then: 

i. For any 02 > 0 , 

2 l;(r(H) >T)< exp ^ 

for all oj such that g{w) > Ai — 02 and any T such that 02 — oiT < 0 . 

ii. Now consider the special case when B is a level set, i.e. suppose for some 04 > 2A2, 
B = {uj : g{w) > Ai — 04} and B' = {u : g{w) < Ai — 204} . Then for all iv G B and 
allT>^ 

ai 

eUt{b') >r')>i- 

wherer = 

The proof of the above lemma follows easily from the standard Azuma-Hoeffding in¬ 
equality for submartingales. We defer it to Appendix C. However we immediately see some 
applications. Recall the dehnition of fig from Table 1 . 

Lemma 11 . Given any e > 0 there exist positive constants c = c{e),d = d{e),N = N{e) 
such that for all n > N and a G Ll 

P^(r(a) > dn^) < e-^^\ ( 71 ) 

Proof. The proof follows from Lemma 10 i. The stochastic process we consider is the com¬ 
petitive erosion chain at and Xt = h{at) is the height function dehned in ( 14 ). We make 


exp 


ai 


32AlT 


-I- exp - 


'siAiTj 


afT^ 


( 68 ) 

( 69 ) 
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the following choice of parameters: 


B = ne 

Ai = 

A 2 = 1 

ai = a(e) appearing in Theorem 5 
02 = 


T = 


2n^ 

Ol 


Clearly (20) and (21) are satisfied by (15) and (16). Thus by Lemma 10 i. 


. 2n^ \ ain- 

T{fle) > - ) < e 

Ol 


□ 


5.1 Proof of Lemma 2 

The proof follows from Lemma 11 and the containment (Remark 1) Pe C T^. □ 


5.2 Proof of Lemma 3 

The proof follows from Lemma 10 ii. Let 

X{t) = h{at) 

where /i(-) is the height function defined in (14). We make the following choice of parameters: 

B = r, 

B' = Ts, 

Ai = 0(1 — aj 2 )n^ 

04 = en^ 

A 2 = 1 

Ol = a(e) appearing in Theorem 5 
T = r?. 


The second containment in Remark 1 along with Theorem 5 satisfy the drift condition (70) 

— ^ 32 ^ — ) • Thus by 

Lemma 10 ii. for all u G Tg, 


(T(r 2 e) > T') > 1 - [exp ( - 


"3^ 


+ exp - 


afn^ 


Thus the proof is complete. 


□ 
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6 Proof of the main result 


The proof of Theorem 1 has the same structure as the proof of Theorem 4. We first state 
and prove the following two results analogous to Lemmas 2 and 3. 

Lemma 12. Given any e > 0 there exist positive constants c = c{e),d = d{e), such that for 
all large enough n and a € 


^a{T{A^) > dn^) < e-"". (72) 

Lemma 13. Given e > 0 there exist positive constants 6,b,d > 0 such that for all large 
enough n and a G n r ^/2 

> e^^) >l-e-^. 

For the dehnition of the sets appearing in the above statements see Table 1. The proofs 
of the above lemmas are significantly more delicate than the proofs of Lemmas 2 and 3. We 
now proceed to the proofs. However for that we need some estimates about internal DLA 
on cylinder graphs. 

6.1 IDLA on the cylinder 

Recall the definition of internal DLA from Subsection 1.1. In this section we discuss the 
case when the underlying graph is the inhnite cylinder Cn x Z>o and the starting locations 
of the particles are uniformly distributed on Cn x {0}. By way of comparison with the 
results stated in [ 8 ], we will show a cruder bound on the fluctuation but we will show that 
the failure probability is exponentially small. We will also need to consider slightly more 
general starting configurations, as described below. 

We now formally define a generalized IDLA process. Consider the graph 

Cn = Cn X Z>o- (73) 

For any integer j > 0 we also dehne the set 

Yj,n = CnX {j}. (74) 

Definition 9. For integers t > 0 the sequence of random subsets (cluster) I{t) is dehned 
inductively. Let 1(0) be any arbitrary subset of Cn- Now given I{t — 1) start a random walk 
uniformly on the set Yq 

I{t) \ I{t - 1) 

consists of the site at which the random walk exits I{t — 1 ) for the first time. 

Remark 4. For the purposes of this article we will consider the initial cluster to be a union 
of rows, i.e. 

7(0) = Cn X In, 

with In = {ii,i 2 , ■ ■ ■ im} for non-negative integers {ii, Z 2 ,... im} with m allowed to depend 
on n. 
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m m m 

Figure 10: The IDLA cluster I{t), and the cluster A{t) = \ /(O)). 


For convenience we introduce the following notation. Let 

(j) : Z>o —>■ Z>o \ Tn 

be the bijection that preserves order. Abusing notation slightly we write 

(j) ■ Cn ^ Cn\ I {0) 

to denote the map that is identity in the first co-ordinate and </> in the second coordinate. 
Also we define 

A{t) ■.= r\m\m)- 

with A(0) = 0 i.e. A{t) is the set of new sites in the growth cluster under the map See 
Figure 10. We now state the result about the growth cluster in the generalized version of 
IDLA on the cylinder to be proved in Appendix A. 

Theorem 6. Given positive number k < 1 for any small enough e there exists a positive 
number c = c{k, e) such that for all large enough n and all In with m < n 

X [0, (1 — €)A:n] C A{kn^) C Cn x [0, (1 + e)kn]j > 1 — 

6.2 Proof of Lemma 12 

We start by defining some notations and making a few remarks. Recall from (10) the 
definition of Bi{a) for i = 1,2. Given cr G D U D' for i = 1,2, let Ii{t) be the blue and 
red IDLA clusters respectively with initial clusters /i(0) = Bi{a). Note that for the process 
l 2 {t) the particles start from x {1}. 

Also for i = 1, 2 let Ii{t) be variants of the processes Ii{t) where the blue and red random 
walks are killed on hitting the lines y = a — ^ and y = a + ^ respectively. We will refer 
to these as the killed IDLA processes. One immediately notices that starting from the same 
initial condition the trivial coupling (using the same random walks for both the processes) 
gives the inclusions 

Ii{t) c 
Bi{at) I 
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hit) 

hit) 


(75) 

(76) 









for i = 1,2 and all times t. Now suppose that and are two initial conditions with 
C i3i(cr^^^). Then again under the trivial coupling 

C (78) 

where for j = 1, 2, and are the blue IDLA and killed IDLA clusters with initial 

cluster Clearly similar statements hold for the red cluster. We omit them for 

brevity. 

Throughout the rest of the proof we will think of the competitive erosion processes 
Bi{at),B 2 {(Tt) and the IDLA processes Ii{-),l 2 {-) starting from Bi{ao), B 2 {cro) respectively 
to be all coupled under the trivial coupling. For neatness we will say that an event holds 
with overwhelming probability if the probability of occurrence of the event is at least 1 — e~^'^ 
for some c > 0 not depending on n. Also let 


T = 

10 


where e appears in the statement of Lemma 12. 


Proof of Lemma 12. By symmetry it suffices to show that at time T with overwhelming 
probability all vertices below the line y < a — y/e are blue. 

Define the event 


Be := 


does not intersect the line y = a 


—, for t < T 
2 


Now given ao = a € 


Be occurs with overwhelming probability. (79) 

To see this notice that there are at least (^ — e)n lines between the levels y = a — ^ and 
y = a that have no red vertices. Thus by the upper bound in Theorem 6 with overwhelming 
probability l 2 {t) will not intersect y = a — ^ for t <T. (79) now follows from (76). 

We now show that on the event Be with overwhelming probability at time T all the initial 
red vertices below the line y = a — in the competitive erosion process become blue, i.e. 

Cn X [0, a - \/e] C Bi(ctt). 

The strategy to show this is to compare the processes Bi{at) and /i(t). Since on the event 
Be no new red particles hit the line y = a — ^ up to time T the blue competitive erosion 
cluster dominates the blue killed IDLA cluster. That is, 

Bi := |7i(t) C Bi{at) for f < t| occurs with overwhelming probability. (80) 

Note however that the domination works in the reverse direction for the actual blue IDLA 
cluster by (76) and hence the need for the killed process. We would be done at this point 
if the lower bound in Theorem 6 was stated for Ii{t). However since Theorem 6 is for Ii{t) 
we need to argue a little more. 
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Now given the initial condition a consider a modification a' obtained by making all the 
rows between the lines y = a — \fe and line y = a — ^ completely red. Then by (77) for 
all t 

Iiit)chit) (81) 

where the LHS and the RHS are the killed blue IDLA clusters with initial conditions a' and 
a respectively. 

Also define /((t) to be the blue IDLA cluster with starting condition a'. By the upper 
bound in Theorem 6 with overwhelming probability /((t) does not intersect the line y = 
a — ^ for t < T. This is because in a' all the vertices between y = a — yT and line 
y = a — ^ are colored red. Note that this implies that both the blue IDLA and the killed 
IDLA clusters starting from a' are exactly the same for the first T steps. That is, 

B 2 := = I[{t) for t < r| occurs with overwhelming probability. (82) 

Now by the lower bound in Theorem 6 

B3:={CnX[0,a-V~e]cl[{T)} occurs with overwhelming probability (83) 

since a' has at most en lines containing red vertices below the line a — y/e. Lastly notice 
that (81) implies the following containment of events: 

n B2 n B3 C jCn X [0, a — -y/i] C Bi(iTr)} • 

Thus we are done since by (80),(82) and (83) all the event on the LHS occur with over¬ 
whelming probability. 

□ 

Having proved Lemma 12 we state and prove a few preliminary lemmas required for the 
proof of Lemma 13. 

Lemma 14. For small enough e > 0 for any cr G n 

P<,(t(D\. 4^) > en2) > 1-e-"" (84) 

for some c = c(e) > 0. For the definitions of Ae and Dg refer to Table 1. 

Proof. As cj G Dg there are at least {^/e — e)n rows with no red vertex between the lines 
y = a — y/e and a and similarly at least (\/e — e)n rows with no blue vertex between the 
lines y = a and a + ^/e. The proof now follows from the upper bound in Theorem 6 which 
implies that the IDLA processes Ii(t),/ 2 (t) require at least en^ rounds with probability at 
least 1 —for some positive c = c(e) before a blue vertex is seen above the line y = a + y/e 
and a red vertex is seen below the line y = a — \/e or a . Thus we are done by (76) which 
says that the same holds for the competitive erosion clusters Hj(t) for i = 1,2. □ 

For any set A <zTl define the positive return time t~^{A) to be 

inf{t > 1 : Ut G A}. 

Lemma 15. Given positive numbers e, c there exist positive constants 6, D such that for 
any u G n T^ 

Po-('r^(De) > cn^) < e~^^^. (85) 
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For the definition of the sets see Table 1. 


Proof. The proof follows from Lemma 10 i. First let a' = ai for the competitive erosion 
chain starting from ag = cr. If a' G we are done. Otherwise by (16) and the fact that by 
hypothesis cr G F^, 

X := h{a) > a{l — al2)v? — 5v?. 

Now let 

X{t) = h{at+i) 

where h{-) is the height fnnction defined in (14). Thus X(0) = x. We make the following 
choice of parameters: 

Ai = a(l — af2)n‘^ 

A2 = 1 

ai = a(e) appearing in Theorem 5 
02 = with 6 < ac /2 

rj-\ _ 2 

1 = cn . 

The choice of 02 works since a G F^ by hypothesis. Also the choice of 6 ensures that 
02 — oiT < 0 and hence the hypothesis of Lemma 10 i. is satisfied. Thus for 5 = ^ by 
Lemma 10 i. for all cr G H F 5 , 

P(^/(r(0£) > T) < exp ^ ■ 

Since starting from cr is one more than r(ne) starting from a' we are done. □ 

Lemma 16. Given e > 0 there exist positive constants 6, b such that for large enough n and 
any oGrienF^nA 

lP<.(r+(a) < r(L! \ A^)) > 1 - 
Proof. Since cr G lie n by Lemma 14 

P,(r(L!\A^)>en 2 )>l-e-'*" 

for some h = h{e) > 0. Now using c = e/2 in Lemma 15 we get that we can choose 8 such 
that for any crGrienF^nA ^ 

Pa(r+(a) > enV2) < (86) 

for some positive constant D. Thus for such a 8 

P^(r+(nj < r(n\ A^)) > P(^(T+(ne) < e/2n^)-IPa(r(n\ A^) < en^) 

> 1-e-^” 

Hence we are done by choosing 6 = |. □ 

We are now ready to prove Lemma 13. 
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6.3 Proof of Lemma 13 


We will specify 5 later. However notice that for any small enough <5 if cr G r ^/2 then by 
Lemma 3 there exist positive constants c, d depending on 6 such that for large enough n 

P,(r' > > 1 - (87) 

where r' = t(H \ F^). Now notice that by hypothesis a Q Ae and hence a G He. Let 
r(i), r( 2 ) ... be successive return times to i.e. r(i) = 0 and for all i > 0 , 

T(i+i) = inf{t : t > G He}- 

The following containment is true for any positive b': Let 


Then 


{r(e4^) < sj C {3i < s such that cr, 


^u}u(y {{. 

\j=o 


Hi) < 


<T-o+i)}n{cT,y, gT^} 


The above follows by first observing whether there exists any i < such that Ui ^ F^. 
Also let j be the first index such that 


Hj) < ^ ^(j+i)- 

Notice that on the event {t(A^) < sj, 

j < s. 

Also by definition on the event {< 7 * G F 5 for all i < s}, 

G F^ n Hf n 

Thus by the union bound, 

P(r(A^) < "■) < P(3 i < s such that ai ^ F 5 ) 

gb'n 

+ er5nH,n.A^}n{r(A-^)<r(,■+!)}). 

i=i 

Recall that by hypothesis a G F^/ 2 - Thus by (87) for any small enough 6 there exists 
d = d{6) > 0 such that for any 6 ' > 0 the first term is less than for large enough n. 

Also notice that by Lemma 16 we can choose 5 such that every term in the sum is at most 
g-fen £qj, gome constant 6 > 0. Thus for such a 6, putting everything together we get that 
for any 6 ' > 0 and large enough n, 

g6 n 

lP(r(A^) < e^'^) < e"'’”. 

i=l 

Hence by choosing h' <h we get that for large enough n, 

P(r(A^) < 


The proof is thus complete. 


□ 


6.4 Proof of Theorem 2 


We will use Lemma 4 and prove something stronger which will be used in the proof of 
Theorem 1. However for brevity we hrst need some notation. We start by recalling the sets 
in Table 1. Now given any positive ei and 62 define the set 

c = d n r 
^€1,62 ^ €- 2 ' 

We claim that for small enough e, (5 > 0 there exists constants D,D' > 0 such that for any 
cr G H, 

K{T{C^,,s/2)>Dn‘^)<e-^'"- (88) 

Clearly this proves Theorem 2 since ^^ , 5/2 d A^. Let 

' t ' = '^(^5/4) 

= r(q/2) 

r'" = r(H,) 

r'"' = r(^^). 

Recall the notational convention made in Subsection 2.2. Now to show ( 88 ) we notice the 
following containment of events: 

[{t < An^} n {t" > e'^"' | 0 ^'} n {t"' < An^ \ Ur'} H {t"" < An^ \ fj,-'"}] C {t(C^ ,5/2) ^ 3Hn^}. 
To see why this containment holds we first notice that 

'^6/A d Vs/2- 

Thus the hrst two events imply that the process hits the set r ^/2 ™ A-n? steps and stays 
inside for an exponential (in n) amount of time, and in particular stays in Vs /2 from time 
Av? through time 3An^. In addition, the third and fourth events together imply that 
regardless of where the chain is at time An?, the chain enters Hg by time 2An^ and then 
enters A^ by time SAn^. Hence in the intersection of the four events, the hitting time of 
C^ s /2 — A^ n Vs /2 is at most 3Hn^. Let D = 3A. Now for a large enough constant A, 
there exists D' > 0 such that the probabilities of of all the four events on the left hand side 
are at least 1 — This follows by Lemmas 2, 3, 11 and 12 respectively. Hence by the 

union bound ( 88 ) follows for a slightly smaller value of D'. □ 

6.5 Proof of Theorem 1 

Note that ( 88 ) is true for all small enough choices of e and 6. However given e by Lemma 
13 there exist 6,b,d > 0 such that for a G C^^s /2 

> e-"") > 1 - e-'". (89) 

For such a 5 the proof of Theorem 1 follows by using Lemma 4 with the following choices 
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of parameters: 


B = -A^i/a 
ti = Dn^ 


t2 

Pi 

P2 


edn 

e-D'n 


The above choice of parameters satisfy the hypotheses of Lemma 4 by (88) and (89). □ 


6.6 Proof of Theorem 3 

The proof is a simple corollary of Theorem 1. Let k = ^. For i = 1.. .k, let 


e 



Now define 

Ai = {{x,y) G Cn : A(x,y) < atn^}. 

Choose 6 = By Theorem 1 for all i = 1, 2,... A: there exists /3i such that 


■kn (V(x, y) G Ai : y < ai +6) > I - e 
kn iy{x,y) ^ Ai : y > Ui - 6) > 1 - 


(90) 


By union bound there exists /3 such that the events on the LHS of (90) simultaneously hold 
for all i with failure probability at most Clearly this completes the proof. To see this 

suppose that all the events hold and there exists (x, y) G Cn such that — y > e. Let 

1 < J < ^ be such that 

aj-in^ < \{x,y) < ajv? 

(take oo = 0). Now by the second event in (90) y > oij-i — 6 which implies that 


^ , X _ ^ , 100 

— y "A oij — oij—\ +0 — o T € 


< e 




which is a contradiction. Similarly one can show 


K^,y) 

ri.2 


y > —e. Details are omitted. □. 


Appendix 

A IDLA on the cylinder 

The proof of Theorem 6 follows by adapting the ideas of the proof appearing in [11]. The 
proof in [11] follows from a series of lemmas which we now state in our setting. Recall (74). 
Let Tz, fkn be the hitting times of and (j){Ykn,n) respectively. 
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Lemma 17. For any z = (x, y) G Cn with y < kn 


kn^^Yo,„iTz < fkn) > ^ ^<f,{w)iTz < fkn) 

w^CnX [ 0 ,kn) 

where 1 P 0 (^) and lPyo,n random walk measures on Cn with starting point 4>{w) and 

uniform over lo,n respectively. 


Proof. By symmetry in the first coordinate under ^ for any j the distribution of the 
random walk when it hits the set Yj^n is uniform over the set Hence by the Markov 
property the chance that random walk hits (p{z) before Y^(^j^n),n after reaching the line 
is 


^ ^ IP'ui(7z ^ ^kn) 


n 


Thus clearly for any j < kn 


^Yonip < fkn) > - ^wiTz < fkn)- 

’ n 

,ri 

The lemma follows by summing over j from 0 through kn — 1. 


(91) 

□ 


Lemma 18. Given positive numbers k and e with e < 1, then there exists (5 = fd{k,€) such 
that for all z = {x, y) with y < {I — e)kn 

^Yoni'^z^ ^ fkn,n) ^ • 

Proof. Since (p{kn) — 4>{y) > ekn the semi-disc of radius min(eA:, l/2)n around z lies be¬ 
low the line Yfj,(^kn),n- The random walk starting uniformly on Yo^n hits the interval {z — 
min(e/c/2, l/4)n, z+miB.{ek/2, l/4)n) with probability at least min(eA:, 1/2). Now the lemma 
follows by the standard result that the random walk starting within radius n /2 has 
chance of returning to the origin before exiting the ball of radius n in Z^. This fact can be 
found in [10, Prop 1.6.7]. □ 

The next result is the standard Azuma-Hoeffding inequality stated for sums of indicator 
variables. 


Lemma 19. For any positive integer n if Xi i = 1, 2 ... n are independent indicator vari¬ 
ables then 

"■ 2 
P(| Xi — p\ > t) < 2e“4;i 

7=1 

n 

where ^ = E Xi. 

2=1 
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A.l Hitting estimates 

Consider the simple random walk {Y{t))t>o on Z^. 

Lemma 20. For {x,y) G 7?, let h{x,y) = x {0})) = (0,0)} be the probability 

of first hitting the x-axis at the origin. Then 


y 1 

) y) 7 2 i ^ 2 I 2 ^■ 

7r(a;"' + y^) x^ + y^ 


Proof. Let 


h{x,y) = < 


Hx,y), y>o 

0 , y = 0 

^-h{x,y), y<0. 


The discrete Laplacian 


Ah{x,y) = h{x,y) - 


h{x + 1, y) + h{x - l,y) + h{x, y + 1) + h{x, y - 1) 


vanishes except when (x, y) 
follows that 


= (0,±1), and A/i(0,±l) = ±\. 
~ a{x,y + l)-a{x,y-l) 

Hx,y) = -^- 


Since h vanishes at oo it 
(92) 


where 


a{x, y) = - log(x^ + y‘^) + K + 0{ \ ) 

TT x^ + y^ 

is the recurrent potential kernel for 7? (see [9] ). Here k is a constant whose value is irrelevant 
because it cancels in the difference (92). □ 

Let A(-) be the simple symmetric random walk on the half-infinite cylinder Cn = Cn xZ>o. 

Lemma 21. For any positive integers j < k, with A = k — j < n: 


i. for any w £ Yfc,n; 

^w{rij) < T+(/c)) < ^ 

where T{j) and T~^{k) are the hitting and positive hitting times ofYj^n and Y^^n respectively 
for X{-). 

ii. there exists a constant J such that for w £ Yj^n o.nd any subset B C Tfc,n; 

P^(A(r(fc)) GB)< J\B\/A. 

Proof, i. is the following standard result about one-dimensional random walk: starting from 
1 the probability of hitting A before 0 is ^. 

Now we prove ii. Clearly it suffices to prove it in the case when B consists of a single 
element. It is easy to check that if Y{t) = (Yfif),Y 2 {t)) is the simple random walk on Z^, 
then 

X{t) = {Yfit)rnodn,\Y 2 {t)\) (93) 
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is distributed as the simple random walk on Cn- For any £ G Z let ti(£) be the hitting time 
of the line y = £, for Y{t). Clearly by (93) for w = (0, j), z = {zi, k) G Cn, 


Fnj{X{T{k)) = z) = '^ P(0j) {Y{Ti{k) A Ti{-k)) G {{zi + in,k), {zi + in, -A;)}} . 

2 = —OO 

By union bound the RHS is at most 

OO OO 

P{o,j){Yin{k)) = izi + in,k)} + ^ P(^o,j) {Y{n{-k)) = (zi + in,-k)} . 

i=—oo i=—oo 

Using the notation in Lemma 20 we can write the above sum as 

OO 

Y. + in,k + j)] 

i=—oo 

By Lemma 20 the above sum is 

Hence we are done. 


□ 


A.2 Proof of Theorem 6 

Equipped with the results in the previous subsection the proof of Theorem 6 will now be 
completed by following the steps in [11] . 


Lower bound: It suffices to show Cn x [0, (1 — e)A:n] C A{{1 + e)kn‘^). Fix G Cn x 
(1 — e)kn For any positive integer i we associate the following stopping times to the 
walker: 


• (T*: the stopping time in the IDLA process 

• r®: the hitting time of (f){z) 

• ''“te.n- hitting time of the set cl){Ykn,n)- 
Now we define the random variables 


N 

M 

L 


{l+e)k'n? 

"Yj ^('r*<crL> the number of particles that visit (t){z) before stopping 
i=l 

(l+tjkrC 

the number of particles that visit (j){z) before reaching (t>{Yknn) 

i=l 

(l+e)/cn^ 

E lt i the number of particles that visit (j){z) before reaching (i>iXknn) 

i=l 

but after stopping. 


Thus 


N > M - L. 



Hence 

i H((l + e)fcn2)) = F{N = 0) < P(M < a) + F{L > a). (94) 

where the last inequality holds for any a. Now by definition 

E(M) = (1 + e)kv? < Tkn^n)- (95) 

We now bound the expectation of L. Define the following quantity: let independent random 
walks start from each w G </>(C'„ x [0, kn]) and let 

L = l(r 2 < Tkn,n for the walker starting at w). 

'!iiSi/>(C„x[0,fcn]) 

Clearly L < L. Hence the RHS of (94) can be upper bounded by F{M < a) + F{L > a). 
Now 

IE(-£) = ^ ^wiTz < Tkn,n)- 

wG(f>{CnX [0,fcn]) 

Also by Lemma 17 and (95) 

(i + |)e(L)<E(M). 

Choose o = (1 + e/4)max( ^^^ ,E(L)) where the (3 appears in Lemma 18. Now using 
Lemma 19 we get 

F{L > a) < exp(—dn) 

F{M < a) < exp(—dn) 

for some constant d = d(e, k) > 0. Thus in (94) we get 

F{M < a) +F{L > a) < 2exp(—dn). 

The proof of the lower bound now follows by taking the union bound: 

E(Cn X [0, (1 — e)A:n] C A((l + e)A:n^) < ^ ^ A((l + e)A:n^)^ 

zSCnX [0,(1— 

< 2exp(—dn) 
zSCnX [0,(1—e)A:?i] 

< exp(—ere). 

where the last inequality holds for large enough re when c is smaller than d. 

Upper bound: In [11] the upper bound is proven by showing that the growth of the 
cluster above level (1 + e)A:re is dominated by a multitype branching process. However here 
we slightly modify the proof to take into account that in our situation the initial cluster 
is not empty. We define some notation. Let us denote the particles making it out of level 
4>iYkn,n) hy WI,W 2 ,... and define 

A{j) = A{wj). 

Choose /co = k{l + y/e)n. We define 

ye,n ■= Yko+i,n- (96) 
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Given the above notation let 

Zi{j) ■= Mj) n Yi^n- 

Define 

MJ) ■= HZi{j)). 

Lemma 22. [11, Lemma 7] There exists a universal Ji > 0 such that for all k,€ € (0,1), 
n > N{k,e) and all positive integers j,i 

Mj) < 

\ lyekn 

We include the proof of Lemma 22 for completeness. However first we show how it implies 
the upper bound in Theorem 6. Let us define the event 

F := Cn X (1 — e)kn C A{kn‘^). 



Now let B = B{k) > 0 be a constant to be specified later. Then 

P(H(fcn^) Cn X [0, k{l + B^)n] n F) < P{Zn'{2ken‘^) > 1) 

< Hn'i’^ken^) 


where n' = {kB)y/en — 1 — k^/en. To see why these inequalities are true first note that the 
set Y,^\n is at height less than k{l + B^/en). Hence the cluster at time kn^ should intersect 
Yn',n to grow beyond height k{l + B^/e)n. However on the event F at most 2ekn^ particles 
out of the first kn^ move beyond height kn. Hence the size of the intersection of Yn\n and 
the cluster is at most Zn'{2ken^). Thus we get the first inequality. The second inequality 
follows trivially from the fact that for a non-negative integer-valued random variable the 
expectation is at least as big as the probability of the random variable being positive. Using 
Lemma 22 we get 


//„/(2A;en^) < 


kn I Ji 


kn 


2ken^ 


1 


kn I Ji 


n' y/ekn 
Aken^ 


k{B — l)y/en y/ekn 


4Ji 


{B-l)k 




k{B—l)^/en 


Thus 

i [0, n] X [0, k{l + Byfe)n\ r^F) <kn 
Now by choosing B such that 4Ji < (B — l)k we are done. 


(B—l)^kn 


□ 


Proof of Lemma 22. The rate at which Yi^n grows is at most the rate at which a particle 
exiting height kn reaches the occupied sites in Thus if X{t) is the random walk on 

Cn defined in (73) then for any m 

pLi{m -M) - fJ-eim) < 

< 


sup Py X{Ty ) G i(m) 

y&ku,n L 


J 


kny/e 


(97) 

(98) 
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where the second inequality follows by Lemma 21 ii. Summing over m = 0,1... j we get 


MJ) < 


J 

kriy/e 


m=0 


Iterating the above relation in ^ with fixed j gives us 

/ 1 

a- 

The lemma follows by using the inequality 




□ 


B Green’s function and flows 

We prove Lemma 6. We start by discussing some properties of the ordinary random walk 
on Cyl„ (defined in (8)). For any v G Cyl„ define 

Gniv) = V—visits to the line y = 0 before riCn x {1})1. (99) 

Lemma 23. For any point {x, y) G Cyl„ 

Gn{x,y) = l-y. 


Proof. Consider the lazy symmetric random walk on the interval [0, n] where at 0 the chance 
that it moves to 1 is | and everywhere else the chance that it jumps is By symmetry 
of Cyl„ in the x-coordinate it is clear that for all (x,y) G Cyl„, 4nGn(x,y) is the expected 
number of times that the above one-dimensional random walk starting from ny hits 0 before 
hitting n. The above quantity is easy to compute and is 4n(l — y). □ 

Remark 5. Thus for any u G O U 

h{(x) = ^ Gn{x,y) (100) 


where h{-) is defined in (14). 

We now define the stopped Green’s function. For any A C Cyl„ and v G Cyl„ define 

Ga{v) = —lE^ [# visits to the line y = 0 before . (101) 

4n 

Lemma 24. Given A C Cyl„ such that A n {Gn x {!}) = 0, for all (x, y) in Cyl„ we have 

GA{x,y) = HA{x,y) - y 


where Ha{-) was defined in (27). 
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Proof. Let yt be the height of the walk at time t < Consider the following telescopic 

series: 

OO 

yT(A<^) -yo = ^{yt +1 - 2 /t)i(r(yi'') > t). ( 102 ) 

t=o 

Notice that since An {Cn x {!}) = 0, t < r^c implies yt < 1. We make the following simple 
observation: 

K[{yt+i-yt)liT{A'^)>t)\:Ft] = | > ^) = 0 ( 103 ) 

where J-t is the filtration generated by the random walk up to time t. Taking expectations 
on both sides of ( 102 ), we get 

^ OO 

^ix,y) [2/r(A=)] - y= >t) = Ga{x, y) 

i=0 

and hence we are done. □ 

Remark 6. Note that the above lemma implies for any (x,y) G Cyl„, 

^{x,y)[yT{An\ — y 

since the Green’s function is a non-negative quantity. 

Next we relate the Green’s function to the solution of a variational problem. The results 
are well known and classical even though our setup is slightly different. Hence we choose to 
include the proofs for clarity. As defined in subsection 4.1 let E denote the set of directed 
edges of Cyl„. 

For any function F : Gyl„ —)• M define the gradient VF : F —)■ M by 

'\IF{v,w) = F{w) — F{v) 
and the discrete Laplacian AF : Gyl„ —>• M by 

AF(n) = F(n)-^ J;F(u;). (104) 

W'^V 

Note that the graph Gyl„ is 4-regular. 

Recall the definition of energy from subsection 4.1. The next result is a standard 
summation-by-parts formula. 

Lemma 25. For any function F : Cyl„ —>■ M 

T(VF) = 4 ^ F{v)AF{v). 

The proof follows by definition and expanding the terms. 

For a subset A C Cyl„ recalling the definition of stopped Green’s function let 

/a := VGa- (105) 

Also recall the definition of divergence (31). 
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Lemma 26. For any {x,y) G A 


div{/^)(x,!,) = | » 

Proof. For any v = {x,y) € A by definition 

d[v{fA)iv)=4AGAiv) = iGAiv)-^GA{w) = I ^ ftL7Jise. 

W'^V ^ 

The last equality follows by the definition of Ga in (101) by looking at the first step of the 
random walk started from v. □ 

We now prove that the random walk flow /a on a set A is the flow with minimal energy. 

Lemma 27. 

£ifA) = if£if) 

where the infimum is taken over all flows from {Cn x {0}) f]A to A^ such that for (x, y) £ A 

Proof. The proof follows by standard arguments, see [12, Theorem 9.10]. We sketch the 
main steps. One begins by observing that the flow fA satisfies the cycle law, i.e. the sum 
of the flow along any cycle is 0. To see this notice that for any cycle 

^ 1 ? ^ 2 ? • • • 


where x[s £ Cyl„, 


k—1 k—1 

fAixi, Xi+i) = y^(Gyi(xj+i) - GA{xi)) = 0. 

i=l i=l 

The proof is then completed by first showing that the flow with the minimum energy must 
satisfy the cycle law, followed by showing that there is an unique flow satisfying the given 
divergence conditions and the cycle law. 

□ 


Now suppose A C Cyl„ \ {Gn x {1}). Then 
£{fA) = £{VGa) = Ga{v)4AGa{v) = ^ E ^ ^ 

iiGyl k^Cfi k^Cfi 

The hrst equality is by definition. The second equality follow from Lemma 25 and the fact 
that Ga is 0 outside A. The third equality is by (106). The last equality is by Lemma 24 
since by hypothesis A n {Gn x {!}) = 0. 

Proof of Lemma 6. The proof now follows from (107) and Lemma 27. □ 
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C Proof of Lemma 10 

We first prove i. Looking at the process uj{t) started from a;(0) = a; we see by (70) that the 
process 

= ^Mr(S) - C[t A t{B)] 

with Xq = g{oj) is a submartingale with respect to the filtration Tt- Also by hypothesis 
\Zt+i — Zt\ < 2A2. 

Now by the standard Azuma-Hoeffding inequality for submartingales, for any time t > 0 
such that 02 — ait < 0 we have 

¥ {Zt - Zq <-{ait - a 2 )) < e . (108) 

Let T be as in the hypothesis of the lemma. We observe that the event {Xq > Ai — 02} H 
{t{B) > T} implies that 

Zt — Zd < —{aiT — 02 ). 

This is because by hypothesis Zq = Xq > Ai — 02- Hence on the event t{B) > T 

Zx — Zq = Xx — aiT — Xq < 02 — aiT 
since Xx < Ai by (68). Thus by (108) 

_ (“l'r-a2)^ 

P(r(H) > T I Ao > Ai - 02) < e . 

To prove ii. let loq = u{t{B‘^)). By hypothesis 

X := /(wo) > Ai - 04 - A 2 . 

since by (69) the process cannot jump by more than A2. Clearly it suffices to show 

ai a\T 

^u^^{t{B') >T')>1- + e~^. 

Now consider the submartingale 

Wt = XtAr'Ar" - Ol [t A t' A t"] 

with Wq = X, where t' = t{B) and t" = t{B'). We first claim that 

IPajo {t' ^ >T) < e~ i 6 c 2 t . (109) 

To see this notice that by the Azuma-Hoeffding inequality it follows that 

2^2 

1 

F{Wx-Wo<-aiT/2) (110) 

On the other hand the event t' At" > T implies 

Wx ^ Ai — 04 — aiT 
Wq > A4 — 04 — A2. 
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Thus the event r' /\t" > T implies 


Wt-Wq< 


aiT 


since by hypothesis T > (109) now follows from (110). 

Now on the event {r' A r" < T} P|{t" < r'} , 

Wt < ~ 204 

Wq > Ai — 04 — A2- 

Hence 

{r' A r" < T} n {r" < r'} ^ Wt - Wq <-"j 
since by hypothesis 04 > 2 A 2 . Thus by (110) we have 

“4 

T ({r' A r" < T} ^^{t" < r'}) < e ^. 

Observe that 

T(r' < r") > P(t' A r" < T) - P (^{r' A r" < T} ^{r" < r'}) . 

This along with (109) imply that t” stochastically dominates a geometric variable with 
success probability at most 

al ajT 

Thus we are done. □ 


Future directions and related models. 

Fluctuations. This article establishes that competitive erosion on the cylinder forms a 
macroscopic interface quickly. A natural next step is to find the order of magnitude of its 
fluctuations. Theorem 1 only shows that the flnctnations are o(n). 

Randomly evolving interfaces. Competitive erosion on the cylinder models a random 
interface flnctnating around a fixed line. It can also model a moving interface if the measures 
Hi and H 2 are allowed to depend on time. An interesting example is /ii(t) = Sz^t) where Zi 
and Z 2 are simple random walks independent of everything else in the process: that is, red 
and bine walkers are alternately released from a red sonrce and blue source that themselves 
perform random walk on a slower time scale. 

Another model of a randomly evolving interface arises in the case of fixed but equal 
measures Hi — Figure 11 shows a variant of competitive erosion in the square grid 
Z?. Initially all sites are colored white. Red and blue particles are alternately released 
from the origin. Each particle performs random walk until reaching a site in — {(0, 0)} 
colored differently from itself, and converts that site to its own color. Particles that return 
to the origin before converting a site are killed. One would not necessarily expect any 
interface to emerge from this process, but simulations show surprisingly coherent red and 
blue territories. 
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10 ® 10 ^ 10 ^° 


Figure 11: A variant of competitive erosion on 1? with = 1^2 = <^o- From left to right, 
the set of all converted sites after n = 10®, 10®, 10^® particles have been released. Most red 
particles convert blue sites and vice versa, so that a relatively small number of sites are 
converted. 

Conformal invariance. Our choice of the cylinder graph with uniform sources on the 
top and bottom is designed to make the function g in the level set heuristic (see (5)) as 
simple as possible: g{x,y) = 1 — y. A candidate Lyapunov function for more general graphs 
is 

h{S) = g{v) 

whose maximum over S' C IL of cardinality k is attained by the level set (7). 

A case of particular interest is the following: Let V = D r\ where L> is a bounded 

simply connected planar domain. We take yi = 6z^ for points zi,Z 2 & D adjacent to D^. As 
the edges of our graph we take the usual nearest-neighbor edges of ^7? and delete every 
edge between D and D^. In the case that D is the unit disk with zi = 1 and Z 2 = —1, the 
level lines of g are circular arcs meeting dD at right angles. The location of the interface for 
general D can then be predicted by conformally mapping D to the disk. Extending the key 
Theorem 5 to the above setup is a technical challenge we address in a subsequent paper. 
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